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MASS DIFFUSIVITY EFFECTS ON DROPLETS IN FILM BOILING 


by Glen J. Schoessow* and Kenneth J. Baumeister** 

ABSTRACT 

This paper presents experimental film boiling data for water droplets 
vaporizing into air, argon, nitrogen, helium, and steam. A dimensionless 
grouping is presented which compares vapor removal by diffusion from the 
top of the drop to vapor removal by film boiling from the bottom of the drop. 
This dimensional grouping is in good agreement with the experimental data. 

INTRODUCTION 

Current research in droplet film boiling on hot surfaces is concerned 
with such problems as fuel vaporization in spark- ignition or diesel engines 
(ref. 1); spray cooling or quenching of heated surfaces (refs. 2 and 3); 
liquid droplet removal in mist sections of boilers (ref. 4); and in other 
diverse problems. References 1 to 12 contain a comprehensive current 
literature summary of the recent work in droplet film boiling as well as a 
discussion of the phenomenon. 

Heat-transfer coefficient correlations for drops in film boiling have 
been based mostly on experimental data for drops exposed to air; how- 
ever, in applications as above the drops usually vaporize in a saturated 
atmosphere of their own vapor. The present paper presents experimental 
film boiling data for water drops vaporizing into steam, air, argon, nitro- 
gen and helium at atmospheric pressure. 

Professor Bell (ref. 5) points that at the present time, we simply do 
not know enough about the temperature and velocity fields surrounding the 
drop to accurately assess diffusivity effects. 'If the pure vapor from under- 
neath the drop rises around the drop and enshrouds the upper surface 
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(fig. 1(a)), then diffusive mass transfer can be neglected (though bulk 
evaporation could still take place). On the other hand, if the vapor from 
underneath the drop disperses radially outward away from the surface 
(fig. 1(b)), and the upper surface is in an air atmosphere, then molecular 
and possibly convective diffusion could occur. The fact that the droplet 
upper regions are usually a few degrees below saturation temperature 
tends to support the latter supposition. Wachters (ref. 6) interprets his 
observation that a droplet in air evaporates faster than a droplet in a 
saturated atmosphere, as indicatir^ that diffusive mass transfer exists 
in the former case. ” 

A flow visualization technique developed in reference 15 gives further 
credence to the two- region flow model. If a liquid nitrogen droplet is 
placed on a water surface instead of the conventional solid surface (see 
fig. 2), the liquid nitrogen drop will still enter a conventional film boiling 
state. Heat will be transferred to the liquid nitrogen droplet and the drop 
will begin to evaporate. However, as the nitrogen vapor moves along the 
water surface small amounts of water will diffuse into the nitrogen vapor 
and condense. 

The flow field outlined by the condensed water vapor is quite similar 
to the two-region flow model in figure 1(b). However, the flow patterns 
seen by this visualization technique will not be exactly similar to that 
which occurs in conventional boiling. In conventional boiling, the vapor 
will generally be much hotter. In this case, bouyancy effects will tend to 
reduce the size of region 2. Nevertheless, this flow visualization tech- 
nique is one more additional justification for the two-region model. 

A movie of the developing vapor flow-field is available in a movie 
supplement to reference 15. Figure 3 shows a series of photographs 
taken from the movie. As seen in these photographs, the upper surface 
of the drop is not engulfed by the vapor; thereby leavii^ the upper surface 
of the drop exposed to allow diffusion and convection to occur. 

Baumeister and Schoessow in reference 13 conceded the importance 
of diffusive effects at low plate temperature especially if the droplet is in 
the metastable Leidenfrost state (ref. 14). However, at high plate tem- 
peratures normally associated with droplet film boiling, they argued that 
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diffusive evaporation from the upper surface will be negligible compared to 
boiling evaporation from the lower surface. Their justification was the 
agreement of the experimental data with their generalized theory which ne- 
glected diffusion. Their motivation for neglecting diffusion was the relative 
ease in obtaining a generalized closed form solution for the drop vaporiza- 
tion time. 

If neglecting diffusion is a reasonable assumption for high temperature 
droplet film boiling in air, would the assumption also hold for droplet 
vaporization in nitrogen or another inert gas ? At what temperature level 
in air would neglecting the diffusive effect lead to large errors ? To answer 
these and other similar questions, a dimensionless grouping will be derived 
which compares vapor removal by diffusion from the top of the drop to vapor 
removal by film boiling from the bottom of the drop. 

In particular, this paper will introduce a dimensionless group which 
can be used as a criteria for predicting when diffusion is important. Next, 
drop vaporization data are taken with various atmospheres and compared to 
the dimensionless group. 


SYMBOLS 


A. 


d 

^fb 


Av 


P 


D 


AB 


F 

Gr 


g 

H 


area 

area of drop in which diffusion occurs 

area of drop in which film boiling heat transfer occurs 

area of drop in which radiative heat transfer occurs 

specific heat at constant pressure of vapor 

diffusion coefficient (see table II) 

radiation view factor 

defined by Py(pL " Py)gL /m evaluated at film temperature 
coefficient of gravity 
defined by X/C AT 

ir 
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H* defined by ^*/C AT 

Ir 

film boiling heat transfer cpefficient 
hr radiative heat transfer coefficient 

h* dimensionless heat transfer coefficient (see table I) 

k thermal conductivity of vapor 

L characteristic length defined as [a/g(pj^ - 

M molecular weight 

Ndc dimensionless parameter defined by equation (11) 

P pressure 

Pr Prandtl number C ii/k evaluated at film temperature 

P_ saturation pressure 

PoQ vapor pressure far removed from drop 

R ‘ gas constant 

r^ effective droplet radius 

r* dimensionless radius (see table I) 

Sc modified Schmidt number, ju/D^g(MP /RTg) with ju. evaluated 
at film temperature 

T temperature 

AT temperature difference 

T plate temperature 

T„ saturation temperature 

s 

t time 

V volume of liquid drop 

q 

V* dimensionless volume of liquid drop, V/L^ 

rate of mass loss by diffusion per unit area 
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M 

Py 


a 

a 


rate of mass loss by film boiling per unit area 
rate of mass loss by radiation per unit area 
latent heat of vaporization 


modified latent heat of vaporization, X ( 1 + 


7 S at' 


-3 


20 


viscosity of vapor 
liquid density 
vapor density 
surface tension 
Stefan- Boltzmann constant 


DIMENSIONLESS PARAMETER 


Gottfried, Lee, and Bell (ref. 8) subdivided a mass balance on a 
droplet in film boiling similar to the following: 


( 1 ) 

where represents the rate of mass loss by diffusion per unit area, 
Wjij and Wj, represents the rate of mass lost from the bottom and sides 
of the droplet by film boiling and radiation per unit area, and, of course, 
the right side of equation (1) represents the rate of change of the volume 
of the liquid drop. The mass loss by film boiling and diffusion are as- 
sumed to be mutually exclusive. The film boiling and radiation and com- 
ponent of mass transfer W||^ and W^. are given as 



( 2 ) 


A 


X 


( 3 ) 
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and can be approximated as 



Dab^ 



(4) 


for a liquid at its saturation temperature and pressure. For this case, 
partial pressure of the liquid in the vapor environment is quite small and 
the log mean partial pressure of the inert gas is equal to the system pres- 
sure. 

Substituting equations (2) to (4) into equation (1) gives 


-Ap,Si.= hjbATAjb+h ATA ^ A , (5) 

^ dt ™ ID r r RTgr^ “ 

Since diffusion is expected to be important only at low AT, the 
radiative component of heat transfer is small compared to the film boil- 
ing component. Consequently, for the purpose of simplicity, radiation is 
now neglected. In this case, equation (5) becomes 


-Ap 


dV ^AB^^s 

^ = hf. ATA.^ + A ^ A , 

^ dt ^ RT r ^ 

so 


( 6 ) 


Dividing both sides of equation (6) by the conductive heat transfer term 
gives 


"^L dV 
% ^TAfb dt 


1 + N. 


DC 


(7) 


where 



^AB^Ps 

^fb 


(8) 
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Equation (7) could now be integrated to determine the total vaporization 
time of the drop subject to the restriction that radiation is small. Of 
course, the relationship between the heat transfer coefficient, heat trans- 
fer area, and drop volume must first be determined before the integration 
is performed. 

In the application of equation (7) to small spherical drops, Gottfried 
and Bell (ref. 12) included the diffusion term which necessitated a numer- 
ical solution. On the other hand, Baumeister, Hamill, and Schoessow 
(ref. 9) related the drop size to volume for both spherical and large drops, 
neglected the diffusion term and obtained a simple closed-form solution 
to the problem. 

A cursory examination of Nj^q (eq. (8)) indicates that for large 
AT’s, NpQ in equation (7) might be safely neglected. However, for small 
AT’s, could become the dominant term in the mass balance (eq. (7)). 
To evaluate Nj^(^ for comparison purposes, expressions for h^ and r^ 
as a function of the droplet volume are required. The expressions for 
hji^ can be written as follows : 




(9) 


hfjj = ^ (Gr Pr h*(V) 

L 


( 10 ) 


where the expressions for the dimensionless parameters r* and h* are 
given in table I for various ranges of the dimensionless droplet volume V*. 
The derivations of these expressions can be found in references 9 and 13. 
For large drop volumes, nonspherical drops, equations (8) and (9) repre- 
sent approximations to the spherical case. 

Substituting equations (9) and (10) into equation (8), taking the ratio of 
Ad/A|b equal to 1, introducing the Schmidt number, and rearranging the 
terms gives 


N - — - 

^DC ~ 


H 

Sc 


Pr3 y/4 1 

Gr H* / 


( 11 ) 
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As a first order approximation, = 1 assumes that the upper half 

of the drop loses mass by diffusion while the lower half by film boiling. 
The value of will be determined for water and various inert 

gases. These values will then be compared to the experimental results. 

In particular, for equal to zero there will be no diffusive evapora- 
tion, while for equal to or greater than 1 diffusive evaporation will 

dominate. For example, a water drop evaporating in a steam atmosphere 
would have a value of equal to zero. 

In expressing in the form given by equation (4), the partial pres- 
sure of the liquid in the vapor was assumed to be small in comparison to 
the saturation pressure of the liquid. In the actual situation, however, 
is proportional to the difference between the two partial pressures; 


W , oc p - p 

d s oo 


( 12 ) 


To evaluate (Pg - P^), two idealized models could be used. 

First, as illustrated in figure 1(a), the saturated vapor could rise up 
from beneath the drop, cover the upper surface of the drop, and to a great 
measure restrict the diffusion process from the top of the drop. In this 
case, P^ would approach the value of Pg. On the other hand, as sug- 
gested by Bell (ref. 5), two distinct flow regions could exist about the 
drop, as shown in figure 1(b). A saturated vapor region could form 
around the edges of the drop (region 1) while the surrounding inert gas 
could form a free convection cell about the top of the drop. In this case, 

« Pg (13) 


and the functional form of equation (4) for would be completely 
justified. 

As mentioned in the INTRODUCTION of this paper, the two region 
model seem to represent the actual physical situation. Consequently, 
equation (13) is valid. 
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ENVIRONMENTAL EXPERIMENTS 


To establish the validity of the assumptions that produced the equation 
for Nj^q (eq. (II)), an experiment should be performed in which the diffu- 
sive evaporation is a significant percentage of the total mass evaporated. 
Figure 4 presents the computed (eq. (11)), values of the Nj^q group for 
a water drop in an atmosphere of helium, nitrogen, argon, and air. Ta- 
ble II lists the values of used in the calculations. 

As seen in figure 4, at a plate temperature of about 360° C, the ratio 
of diffusive evaporation from the top of the drop to conductive evaporation 
from the bottom of the drop should be about 10 percent in air for a water 
drop of about 0. 7 cc. In other words (1 + Ng^) equals 1.1. The parame- 
ter (1 + Njjq) is plotted in figure 4 rather than Nj^^, because the vapor- 
ization rate is proportional to (1 + Nj^q), as seen in equation (7). For 
argon and nitrogen atmospheres (1 + Njj^) increases only by 1 or 2 per- 
cent relative to air. Consequently, the evaporation time of a water drop 
in argon or nitrogen will be nearly equal or slightly smaller than the 
evaporation time of the same drop in air. 

On the other hand, as seen in figure 4, the evaporation time of a 
water drop in helium would be much smaller than in air because of the 
increase in mass diffusion removal from the top of the drop as indicated 
by the larger value of (1 + Nj^q). 

To verify the above predictions, water drops were evaporated in air, 
argon, helium, and a nitrogen atmosphere. The experiments were con- 
ducted in a large cylindrical tank 60 cm in diameter and 60 cm long. Dry 
gas was periodically fed into the tank to keep the moisture level down. In 
addition, the first data point taken for a drop in a completely dry environ- 
ment was always repeated at the end of a data set to check that the mois- 
ture level had not inadvertently risen to a high level. 

The results of the experiment are shown in figure 5, which plots total 
vaporization time over drop volume. Qualitatively, the experimental re- 
sults agree with the theoretical predictions. The vaporization times of a 
water drop at a plate temperature of 360° G were predicted to be nearly the 
same for the air, argon, or nitrogen environment. Also, the vaporization 
time in helium is much smaller as predicted. 
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For example, the 25 percent higher value of (1 + for helium is 
mirrored in the data by an approximately 30 percent decrease in the total 
vaporization time for a 0. 7 cc drop. 

In addition to the previous experiments, 1 cc water drops were evapo- 
rated on a stainless steel hot plate at 376® C in air and with a steam cover. 
The average total vaporization time in air was 352 seconds, while in the 
steam atmosphere it was 377 seconds. This 7 percent decrease in the 
vaporization time in air compares favorably to the calculated value of 
1. 08 for (1 + Njpc). 

Equation (7) can also provide an additional check of the dimensionless 
group First, note that the left side of equation (7) is independent of 

the atmosphere. Thus, if the expression of equation (7) for atmosphere 
No. 1 is divided by the expression for atmosphere No. 2, the resulting 
ratio is 



The ratio S for nitrogen (No. 1) to helium (No. 2) was calculated to be 
0. 84 at 0. 7 cc. The measured ratio of the slopes from figure 5 is 0. 91. 
The ratio S for air (No. 1) to helium was calculated to be 0. 82 while the 
measured value is 0. 75, which is fair agreement. 


CONCLUSIONS 


Based on the two- region flow model, the dimensionless group 



Sc 


Gr H*j h|rj 


gives the ratio of diffusive evaporation rate from the top of the drop in 
comparison to film boiling evaporation from the thin vapor gap beneath 
the drop. This group is in fair agreement with the experimental data. 
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TABLE I. - RADIUS AND CONDUCTIVE HEAT TRANSFER COEFFICIENT 

FOR DROPLET FILM BOILING 


Dimensionless 

volume, 

V* = X 


Drop shape 


Dimensionless 

radius, 

r* = L 
° L 


Dimensionless 
heat transfer 
coefficient, 

h? = !l 


^(GrPrH*)^/^ 


V* < 0.8 


Small spheroid 




r* 
0 


*\l/3 


4tt 


h* = 1.1 
s 


0.8 <V* < 155. 


Large drop 




♦ _f l.25 V 

° V 77 


♦5/6>f/2 


) 


h* = 1.075 
s 


V* > 155. 


Extended drop 
(^Constant ^ 
\thickness j 


r* =(MiX 
0 \ 

\ ^ 


. 1/2 


h* = 1.64 
s 
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TABLE n. - DIFFUSION COEFFICIENTS (373 K 
AND ATMOSPHERIC PRESSURE) 


Liquid-gas combination 

^AB’ 

cm^/sec 

Refer- 

ence 

Comments 

Water-air 

0.353 

16 


Water-argon 

0.41 

17 

(a) 

Water- helium 

.. . .. . .j 

1.22 

17 

(b) 

Water-nitrogen 

0.40 

17 

(b) 


^Dab calculated by Slattery- Bird equation corrected 
for water (eq. 11-27) as given in ref. 17. 

•jU 

^Experimental data given in table 11-4 of ref. 17 
was linearly extrapolated to 373 K. 
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Figure L - Velocity flow models for drops in film boiling. 



Figure 2. - Flow visualization technique. 
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Figure 3. - Photographs of vapor flow patterns for nitrogen drop on 
water. 
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Figure 4. - Comparison of (1 + Nq^) for vari- 
ous water and inert gas combinations for 
atmospheric pressure and a plate tempera- 
ture of about 360 C. 



Figure 5. - Total vaporization data for various water and 
inert gas combinations. 
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